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Abstract 



We analyze experimental data for the production of A baryons in e^e~ annihilation in 
terms of scale dependent, QCD evolved, A fragmentation functions. Apart from the vast 
majority of the data for which the polarization of an observed A was not determined, 
we also consider the recent LEP measurements of the longitudinal polarization of A's 
produced on the Z-resonance. Such data correspond to spin-dependent fragmentation 
functions for the A. We point out that the present data are insufficient to satisfactorily 
fix these. We therefore suggest several different sets of fragmentation functions, all com- 
patible with present data, and study the prospects for conceivable future semi-inclusive 
deep-inelastic scattering experiments to discriminate between them. We provide the com- 
plete next-to-leading order QCD framework for all the processes we consider. 



1 Introduction 



Measurements of rates for single-inclusive e'^e annihilation (SIA) into a specific hadron 
H, 

e+e-^{j,Z)^HX , (1) 

play a similarly fundamental role as those of the corresponding crossed "space-like" deep- 
inelastic scattering (DIS) process ep e'X. Their interpretation in terms of scale- 
dependent fragmentation functions D^{z, Q'^), the "time-like" counterparts of the parton 
distribution functions fuiXiQ"^) of a hadron iJ, provides a further important, comple- 
mentary test of perturbative QCD. In analogy with the "space-like" case, Dj{z,Q'^) is 
the probability at a mass scale Q for finding a hadron H carrying a fraction z of the 
parent parton's momentum. QCD completely predicts the Q^-dependence of the process- 
independent fragmentation functions D^{z,Q'^) via the Altarelli-Parisi evolution equa- 
tions, once a suitable non-perturbative hadronic input at some initial reference scale 
has been determined from the data. So far, only the fragmentation into the most copiously 
produced light mesons (vr, K) has been the issue of a thorough QCD analysis 

It is one of the purposes of this work to study along similar lines as in [Q whether 
such a formalism also applies to the production of A-baryons. A's are also produced 
at fairly large rates and, as for pions and kaons, the range covered by present SIA 
experiments p|, |^ allows a detailed quantitative QCD analysis. Recently measured A 
production rates in semi-inclusive DIS (SIDIS) provide an important testing ground 
for the fragmentation functions extracted from SIA data. 

The production of A baryons appears to be particularly interesting also from a dif- 
ferent point of view. Contrary to spinless mesons like pions and kaons, the A baryon 
offers the rather unique possibility to study for the first time spin transfer reactions. The 
self- analyzing properties of its dominant weak decay A p7c~^ and the particularly large 
asymmetry of the angular distribution of the decay proton in the A rest-frame allow 
an experimental reconstruction of the A spin. Over the past years "spin physics" has 
attracted an ever growing interest, as experimental findings have not always matched 
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with "naive" theoretical expectations, the Gourdin-Elhs-Jaffe sum rule [0 being the most 
prominent example here. Studies of A polarization could provide a completely new in- 
sight into the field of spin physics whose theoretical understanding is still far from being 
complete despite recent progress, and they might also yield further information on the 
hadronization mechanism. 

In a strategy was proposed for extracting in SIA the functions ADj (2:, Q^) describ- 
ing the fragmentation of a longitudinally polarized parton into a longitudinally polarized 
A i, 

ADf{z,Q')^D^l^]{z,Q')-D^[;]iz,Q') , (2) 

where D^^'^^{z,Q'^) {D^j-^^ {z , Q'^)) is the probability for finding a A with positive (neg- 
ative) helicity in a parton / with positive helicity (by taking the sum instead of the 
difference in (^ one recovers the unpolarized fragmentation function -Dj). If the energy 
is far below the Z-resonance, one longitudinally polarized beam is required in order to 
create a non- vanishing net polarization of the outgoing (anti) quark that fragments into 
the A, and to obtain a non-zero twist-two spin asymmetry. At higher energies, such as at 
LEP, even no beam polarization is required since the parity-violating qqZ coupling auto- 
matically generates a net polarization of the quarks. Here, ALEPH [|^, DELPHI |[TT| , 



and OPAL |[T2[ have recently reported first results for the polarization of A's produced on 



the Z-resonance. 

Realistic models for the ADj^{z,Q^) are also of particular relevance for reliable esti- 
mates of production rates and spin transfer asymmetries at present and future dedicated 
spin experiments. Here the AD^^ZjQ"^) can be probed in SIDIS or photoproduction in 
the current fragmentation region, Ip —>■ I' AX, where either a longitudinally polarized lep- 
ton beam or a polarized nucleon target would be required. Such measurements can be 
carried out at HERMES |13[ and are planned by the COMPASS collaboration. After 



the scheduled upgrade of the HERA electron ring with spin rotators in front of the HI 
and ZEUS experiments, longitudinally polarized electrons will be also available for high- 
energy ep collisions, and similar measurements with polarized A's in the final state could 
be performed here. Furthermore, having also a polarized proton beam available at HERA 
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p!5| would allow the measurement of various different twist-2 asymmetries, depending on 
whether the e and/or the p beam and/or the A are polarized, i.e., ep AX, ep —>■ AX, 
and ep AX (as usual, an arrow denotes a polarized particle). 



So far estimates for future A experiments have relied on simple models ||T6[ or on 
Monte-Carlo simulations tuned with several parameters and parametrizations of scale- 



mdependent spin-transfer coefficients Cf which link longitudinally polarized and unpo- 



larized fragmentation functions via ||T7 



ADjiz)=C}iz)Djiz) . (3) 

Different phenomeno logical models for the exist. A first one is based on the naive 
non-relativistic quark model where only s-quarks can contribute to the fragmentation 
processes that eventually yield a polarized A. Another approach goes back to estimates 
by Burkardt and Jaffe P, for a fictitious DIS structure function of the A, for which 
sizeable negative contributions from u and d quarks are predicted by analogy with the 
breaking of the Gourdin-Ellis- Jaffe sum rule ^ for the proton's Qi- It is then assumed 
that such features also carry over to the "time-like" case [1^ (see also ^7^). Of course, 
relations like (^ cannot in general hold true in QCD. Due to the different Q^-evolutions 
of and Dj, it cannot be correct to assume scale independence of the Cj in (|^), and 
therefore one has to specify a scale at which one implements such an ansatz. 

It is the main purpose of this paper to address the issue of fragmentation into polarized 
A's in a detailed QCD analysis. Here it will be possible for us to work even at next- 
to-leading order (NLO) accuracy, as the required spin-dependent "time-like" two-loop 
evolution kernels were derived recently For the first time, we will provide some 

realistic sets of unpolarized and polarized fragmentation functions for A baryons. A 
useful restrictive constraint when constructing models for the is provided by the 

positivity condition (similarly to the "space-like" case), i.e. 



\ADj{z,Q')\<Djiz,Q') , (4) 

with the D'j{z,Q'^) taken from the unpolarized analysis. As the available sparse data 
from LEP are by far not sufficient to completely fix the ADj^(z,Q^), we will propose 
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several different sets of polarized fragmentation functions, all compatible with the LEP 
data. Some of the sets will be based on the ideas outlined in the previous paragraph. 
Our various proposed are particularly suited for estimating the physics potential 

of future experiments to determine the polarized fragmentation functions more precisely. 
We hence present detailed predictions for future SIDIS measurements at HERMES and 
the HERA collider. In this context we also provide the necessary framework to calculate 
helicity transfer cross sections in SIDIS at NLO. 

The remainder of the paper is organized as follows: in the next Section we develop 
the formalism for unpolarized SIA and discuss in detail our analysis of leading order (LO) 
and NLO A fragmentation functions. In Section 3 we turn to the case of longitudinally 
polarized A production and present our different conceivable scenarios for the ADf{z,Q^). 
In Section 4 we compare our unpolarized distributions with recent A production data in 
SIDIS and study the potential of present (HERMES) and future (HERA) spin physics 
experiments to discriminate between the different proposed sets of polarized fragmentation 
functions. Finally our results are summarized in Section 5. The Appendices collect the 
required unpolarized and polarized NLO coefficient functions for SIA and SIDIS. 



2 Unpolarized A Fragmentation Functions 



In the last few years several experiments Q have reported measurements of the unpo- 
larized cross section for the production of A baryons, which allows an extraction of the 
unpolarized A fragmentation functions required for constructing the polarization asym- 
metries and as reference distributions in the positivity constraint (^). We emphasize at 
this point that the wide range of c.m.s. energies covered by the data |0, ^ (14 < ^/s < 91.2 
GeV) makes a detailed QCD analysis that includes the Q^-evolution of the fragmentation 
functions mandatory. 

The cross section for the inclusive production of a hadron H with energy Eh in SIA at 
a c.m.s. energy ^/s, integrated over the production angle, can be written in the following 
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way pg, pTI: 



1 da 



H 



(Ttot dxE el 



[2Fl'{xE,Q')+Fl!{xE,Q' 



(5) 



where xe = '^Ph ■ i/Q"^ = '^Eh/^/s (g being the momentum of the intermediate 7 or Z 
boson, = = s) and 



CTtot 



47^a2(Q2^ 



1 + 



asm 



7T 



(6) 



is the total cross section for e~^e~ —>■ hadrons including its NLO 0{as) correction. The 
sums in run over the nj active quark flavours g, and the Cq are the corresponding 

appropriate electroweak charges (see Appendix A for details). 

To NLO accuracy the unpolarized "time-like" structure functions and F^ in (^) 
are given by 



2Ff{xE,Q') = 5^e;|[Df(a;^,Q2) + Df(x^,Q^ 



(7) 
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[Cl ® {Df + Df) + C] ® ] [xE. Q' 



FE{xe.Q') = ^E^? [Ci^iD^ + Df) + C^0Df]ixE,Q') , (8) 



27r 

q 

with the convolutions ® being defined as usual by 

{C®D){XE,Q')= ^c(^)D{y,Q^^ 



(9) 



The relevant NLO coefficient functions C^'g ^ can also be found in Appendix A. 

To determine the Q^-evolution of the in Eqs. (P) it is as usual convenient 
to decompose them into flavor singlet and non-singlet pieces by introducing the densities 



D^. and the vector 



where 



H 



D 



H 



D 



DE^j:iDf + Df) 



(10) 
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One then has the following non-singlet evolution equations {q, q being two different fla- 
vors) : 



d 



cilnQ2 



H 



qq 

DiT) 

qq,- 



D 



H 



(12) 
(13) 



The two evolution kernels Pqq'±{z, as{Q^)) become different beyond LO as a result of the 
presence of transitions between quarks and antiquarks. The singlet evolution equation 
reads 

d 



P(^) ® D 



H 



d\nQ^' 

where we write the singlet evolution matrix P^'^^ as: 



Pc 



(T) 



qq 



2nfP^P 



1 piT) 
2nf ^19 



Pc 



(T) 



(14) 



(15) 
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To NLO, all splitting functions [22 



PiP{z,a,)={ 



in (p!2D-(p!5|) have the perturbative expansion 

(z) . (16) 



It should be noted that the evolution equations can be straightforwardly solved analyti- 
cally in Mellin-n space along the lines as, e.g., described in |21|. The desired Dj{z,Q'^) 
are then obtained by a standard numerical Mellin-inversion. Needless to mention that the 
corresponding LO expressions are entailed in Eqs. (§)-(§) by simply dropping all 0{as) 
contributions and by evolving the P'f{z, Q"^) in LO. 

In our numerical analysis we cannot include data with xe < 0.1, for two reasons: first 
of all, for small xe-, finite-mass corrections to Eq. (|5D proportional to AMj^/ sx\ become 
more and more important, but are not accounted for in the calculation. There is also a 
more severe limitation set by the evolution equations outlined above: the NLO "time-like" 



splitting functions Pgq 



z) and P( 
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in (|I5]) turn out be much more singular than 
their corresponding "space-like" counterparts as z — > 0. While the leading small-x terms 
in the "space-like" case are proportional to l/x, 

(17) 



lhnP(f(x) 
limP(f(x) 



as 
2^ 



2Ce as 9CeCa - 40C^T. 
h 



X 
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9x 



as (2Ca ^ as UCeTj - AQCaTj 
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the "time-like" splitting functions show an even stronger negative behaviour for z ^ 
due to the dominant large logarithmic piece ~ In^ z/^; in the NLO part, 



(see for example |^), with the usual QCD colour factors Ca = 3, Cp = 4/3, and 
Tf = T^Uf = nf/2. This singular behaviour of the "time-like" splitting functions is even 
so strong that it may ultimately lead to negative NLO fragmentation functions in the 
course of the Q^-evolution and hence to unacceptable negative cross sections at some value 
x^; <^ 1, even if the evolution starts with positive distributions at the initial scale. Clearly, 
the description of fragmentation processes by perturbative QCD without resummation of 
small-z logarithms breaks down for values of xe where this happens, and in order to avoid 
these severe problems we include as usual (see, e.g., [|l|) only data with xe > 0.1 in our 
analysis^. 

As pointed out in [Q, the QCD formalism is strictly speaking only applicable to 
strongly produced A's. A certain fraction of the data 0, Q will, however, consist of 
secondary A's resulting from e~^e~ — > with the subsequent decay S° A7, not to 
be included in the fragmentation functions 0. For simplicity, we will ignore this problem 
and (successfully) attempt to describe the full data samples by fragmentation functions 
that are evolved according to the QCD Q^-evolution equations. 

Unless stated otherwise, we will refer to both A° and A°, which are not usually dis- 
tinguished in present e+e~ experiments @, Q, as simply "A". As a result, the obtained 
fragmentation functions always correspond to the sum 



This also considerably simplifies the analysis, since no distinction between "favoured" and 
"unfavoured" distributions is required. Since no precise SIDIS data are available yet, it 
is not possible to obtain individual distributions for all the light flavours separately, and 
^We do not include data either that have been averaged experimentally over a large bin of xe- 




(20) 



(19) 



DfixE, Q') ^ DfixE, Q') + Df\xE, Q') 



(21) 
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hence some sensible assumptions concerning them have to be made. Employing naive 
quark model SUf{3) arguments and neglecting any mass differences between the u, d, 
and s quarlcs, we assum^ that all the light ffavours fragment equally into A, i.e. 

Dt = = Df = D^ = = Df = ■ (22) 
Needless to say that the q and q fragmentation functions in (P^) are equal due to Eq. 



At variance with the usual DIS case, c^zscontinuous heavy quark (HQ) fragmentation 
functions should be included at each heavy flavour threshold (see also Anyway, the 
inclusion of the HQ contributions essentially only leads to a change in the normalization of 
the light quark densities, which would just be larger if the HQ ones were not present. In our 
analysis we start the evolution of the HQ contributions at the mass of the corresponding 
HQ, but the precise value for this is anyhow irrelevant since all the data are in a region 
of s > ml {h = c,b). 

For our analysis, we choose to work in the framework of the "radiative parton model" 
which is characterized by a rather low starting scale /i for the Q^-evolutions. The "ra- 
diative parton model" has proven phenomenologically successful in the "space-like" case 
for both unpolarized and polarized parton densities, and also in the "time-like" 



situation for photon fragmentation functions [^, [27 



At the initial scale (/i|o ~ 0.23 GeV^, /i^^Q = 0.34 GeV^) we choose the following 
simple ansatz: 

Df{z,fi^) = NfZ^'fil- z)'^f , (23) 

where f = q,c,b, g and, as stated, in the case of heavy quarks yU^ = m^. Utilizing Eq. (p2| ) 
and assuming for simplicity that = Nf, = Ng, a total of 10 free parameters remains 
to be fixed from a fit to the available 103 data points p|, Q (after applying the xe cut 
mentioned above). The total values are 103.55 and 104.29 in NLO and LO, respectively, 
and the optimal parameters in (p3D can be found in Table 1. It should be noted that by 



^Fits allowed to be more general do not seem to improve the final /d.o.f. 
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Parameter 


LO 


NLO (MS) 




0.63 


0.55 




0.23 


0.22 




1.83 


2.16 




0.91 


2.23 


ag 


1.36 


1.86 




3.14 


3.48 




-0.41 


-0.35 




5.66 


6.06 




-0.29 


-0.32 


Pb 


5.01 


5.45 



Table 1: Optimal parameters for the unpolarized fragmentation functions in Eq. (^31). 

taking into account an additional 4% normalization uncertainty for the LEP data 0, 
can be further reduced but without any noticeable changes in the distributions. 

A comparison of our LO and NLO results with the data is presented in Fig. 1, where all 
the existing data 0, || have been converted^to the "format" of Eq.(||). One should note 
that the LO and NLO results are almost indistinguishable, demonstrating the perturbative 
stability of the process considered. Furthermore, there is an excellent agreement between 
the predictions of our fits and the data even in the region of "small" xe which has not 
been included in our analysis. 

Fig. 2 shows our LO and NLO fragmentation functions as specified in Eq. ( p3|) and 
Table 1, evolved to = 100 and 10^ GeV^. As can be seen, the heavy quark fragmen- 
tation functions turn out to be comparable to the light quark ones for small z, whereas 
they are suppressed for ^ > 0.3. It should be also noted that our c and b fragmentation 
functions are also in agreement with recent results from SLD ^ for the c/uds and h/uds 

■^The available data sets |^ are presented in terms of three different variables: xe, Xp, and ^. These 
variables are simply related to each other by Xp = j3xE with j3 = -^Z f — rri^^ / Ejj , and ^ = ln(f/a;p). 
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ratios of A production rates in flavour-tagged Z decays. 



In Fig. 3 we show our fragmentation functions for light quarks and gluons as functions 
of for several fixed values of z. One can observe the importance of the QCD evolution 
of the fragmentation functions which we will use for making predictions for SIDIS at 
values much lower than the ones at which the fragmentation functions were extracted. 

Finally, we investigate the contribution of our A (more precisely A + A, see Eq. (^)) 
fragmentation functions D^{z^Q'^) to the momentum sum rule 

f dzzDf{z,Q') = l . (24) 

Eq. (^) expresses the conservation of the momentum of the fragmenting parton / in the 
fragmentation process, i.e. each parton / will fragment with 100% probability into some 
hadron H. Of course the sum rule (|2l) should be dominated, even almost saturated. 



by the fragmentation into the lightest hadrons such as vr and K mesons. Hence the 



contribution to (p^ due to is expected to be rather small. Indeed we find that in LO 



and NLO the contribution of our light quark (gluon) A + A fragmentation functions to 
the momentum sum rule (p^ ) only amounts to about 2 — 3% (1 — 2%). 



3 Polarized Fragmentation Functions 

Having obtained a reliable set of unpolarized fragmentation functions we now turn to the 
polarized case where unfortunately only scarce and far less precise data are available. In 
fact, no data at all have been obtained so far using polarized beams. The only available 
information comes from unpolarized LEP measurements [10-12] profitting from the parity- 
violating electroweak qqZ coupling. 

For such measurements, done at the mass of the Z boson (Z-resonance), the cross 
section for the production of polarized hadrons can be written as [|[ ^ 



3. 

dVtdxE 2 s 



(25) 
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If, as for the quoted experimental results [10-12], the cross section is integrated over the 
production angle 9, the anyway small, charge suppressed contribution from drops out. 
One can then define the asymmetry 

"Fi^ + Ff/2 ^^^^ 
which corresponds to the "A-polarization" observable measured at LEP. The polarized 

structure functions , g^ and gj^ in (|25|) and (|26|) are given in NLO by the following 

expressions: 

gnxE,Q') = J29',{[^K(^E,Q') + ADf{xE,Q')] (27) 



"^^^'^ [ACl ® (ADf + ADf ) + ACl ® ADf] {xe, Q' 



2n 



9^{XE,Q') = J29M^D^i^E,Q')-ADf{xE,Q' 

[AC'^ ® ( ADf - ADf)] {xE, Q') I 
9^{XE,Q') = ^^J29A^Ci®i^Df-AD^)]{xE,Q') 



(2^ 



2n 



2n ^-^L--</-v— 9 — (29) 

with the convolutions as already defined in (|^). The appropriate effective charges gq and 
g'g as well as the required spin-dependent MS coefficients AC^'I'^ in (^7D -(^9D can be found 
in Appendix B. 

Note that both g^ and gj^ in (|28|) , (^9|) are non-singlet structure functions, and hence 
only the valence part of the polarized fragmentation functions can be obtained from the 
available LEP data [10-12]. In addition the A°'s and A°'s give contributions of opposite 
signs to the measured polarization and thus to g^j^^. Unfortunately, it is clear that the 
available LEP data [10-12], all obtained on the Z-resonance, cannot even sufficiently 
constrain the valence distributions for all the flavours, so some assumptions have to be 
made here. Obviously, even further assumptions are needed for the polarized gluon and 
sea fragmentation functions in order to have a complete set of fragmentation functions 
suitable for predictions for other processes, in particular for SIDIS (see Section 4). 

In the present analysis the heavy flavour contributions to polarized A production are 
neglected, and u and d fragmentation functions are taken to be equal. Furthermore, 
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polarized "unfavoured" distributions, i.e., AD^° = AD^° , etc., and the gluon fragmen- 
tation function AD^ are assumed to be negligible at the initial scale fi, an assumption 
which of course deserves a further scrutiny (we will discuss the impact of choosing a 
different boundary condition for the gluon fragmentation function later). The remain- 
ing spin-dependent quark fragmentation functions are then related to the corresponding 
unpolarized ones taken from Section 2 in the following simple way 

They are subject to the positivity constraints (^), which simply imply a > and |A^„| < 1. 
These input distributions are then evolved to higher via the appropriate Altarelli-Parisi 
equations which are completely similar to the ones presented in Eqs. (p!0|)-(p!6|), with just 
all unpolarized quantities (like, for instance, P^'^^) replaced by their appropriate polarized 
counterparts. For the NLO evolution one has to use for this purpose the spin-dependent 



"time-like" two-loop splitting functions as derived in |T9[ in the MS scheme. Due to the 
rather limited amount of available data it does not appear reasonable for the time being 
to introduce more free parameters than the two in Eq. (^). 

Within this framework we try three different scenarios for the polarized fragmentation 
functions at our low initial scale fi, to cover a rather wide range of plausible models: 

Scenario 1 corresponds to the expectations from the non-relativistic naive quark model 
where only s-quarks can contribute to the fragmentation processes that eventually yield 
a polarized A, even if the A is formed via the decay of a heavier hyperon. We hence have 
N,, = in (|30|) for this case. 



Scenario 2 is based on estimates by Burkardt and Jaffe |T8| for the "space-like" DIS 
structure function of the A, predicting sizeable negative contributions from u and d 
quarks to gi by analogy with the breaking of the Gourdin-Ellis- Jaffe sum rule for the 
proton's g^. Assuming that such features also carry over to the "time-like" case [|18|, we 
simply impose = —0.20 (see also 



Scenario 3: All the polarized fragmentation functions are assumed to be equal here, 
i.e. Nu = 1, contrary to the expectation of the non-relativistic quark model used in 
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seen. 1. This rather "extreme" scenario might be reahstic if, for instance, there are 
sizeable contributions to polarized A production from decays of heavier hyperons who 
have inherited the polarization of originally produced u and d quarks. 



Parameter 


LO 


NLO (MS) 




seen. 1 


seen. 2 


seen. 3 


seen. 1 


seen. 2 


seen. 3 


Nu 





-0.2 


1.0 





-0.2 


1.0 


a 


0.62 


0.27 


1.66 


0.44 


0.13 


1.33 



Table 2: Resulting optimal LO and NLO fit parameters as introduced in for the three 
different scenarios described in the text. 

Our results for the asymmetry in ( p6D within the three different scenarios are 
compared to the available LEP data [10-12] in Fig. 4. The optimal parameters in ( pOD for 
the three models can be found in Table 2. As can be seen, the best agreement with the 
data is obtained within the (naively) most unlikely seen. 3. The differences occur mainly 
in the region of large Xe, where seen. 1 and 2 cannot fully account for the rather large 
observed polarization. It turns out that this is a consequence of the assumed SU{3)f 
symmetry for the unpolarized fragmentation functions, and of the positivity constraints 
(^): for instance, SU{3)f symmetry of the implies in the case of seen. 1 that the 
asymmetry at large x^; behaves asymptotically roughly like —AD^/3D^. Thus, even when 
saturating the positivity constraint (H) at around xe = 0.5 it is not possible to obtain a 
polarization as large as the one required by the ALEPH and OPAL data |T2| . We note 
that the assumed SU (3) / symmetry for the unpolarized fragmentation functions could of 
course be broken. It is clear at this point that further information on the polarized and 
the unpolarized A fragmentation functions is needed, which could be provided by future 
precise SIDIS measurements. 

Finally, in Fig. 5 we show the LO and NLO partonic fragmentation asymmetries for 
each flavour distribution separately, i.e. Aj = /\D^/D^. A positive polarized gluon 
fragmentation function has built up in the Q^-evolution in spite of the vanishing input 
at yU^. In order to analyze the effect of imposing a different boundary condition for 
the polarized gluon fragmentation function, we include in Fig. 5 also the results of a 
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LO fit similar to the one performed within scenario 1, but now using the maximally 
allowed polarized gluon input ADg(yU^) = Dg{ii^) instead of ADgi^fi"^) = 0. Besides the 
expected result of having now a larger gluon polarization at = 10 GeV^, an important 
enhancement of the u and d distributions (which are practically vanishing in the original 
scenario 1) can be observed, which is due to the perturbative generation of sea by polarized 
gluons in the course of the evolution. In fact, at small values of z, the u and d distributions 
become even larger than the ones of scenario 3. Obviously, only different combined further 
measurements, like in e~^e~ annihilation and SIDIS with polarized beams, will be capable 
of determining the gluon (and also the sea) fragmentation function more precisely. 



4 A Production in SIDIS 



Equipped with various sets of polarized fragmentation functions, let us now turn to the 
SIDIS process eA^ — ^ e'HX which should be very well suited to give further information 
on fragmentation functions. In this case, the cross section is proportional to a combination 
of both the parton distributions of the nucleon and the fragmentation functions for the 
hadron H. The latter thus automatically appear in a constellation different from the one 
probed in e'^e~ annihilation. 

In the particular case where both nucleon and hadron are unpolarized, the cross section 
can be written in a way similar to the fully inclusive DIS case pOl, pi], |30|: 



rfa^ 2 7r«2 p(i^(l_^)2) 2{l-y) ^/H, ^ ^2, 



dx dy dzu 



-2 Fi''" (x, ZH, Q') + Z^^_lLF'y" (a;, ZH. Q') 



y y 



(31) 

with X and y denoting the usual DIS scaling variables (Q^ = sxy), and where [ P0| , pT[ | 
Zh = Ph ■ Pn/Pn ■ Q with an obvious notation of the four-momenta, and with — = 
Q^. Strictly speaking, Eq. (|3TD and the variable zh only apply to hadron production in 
the current fragmentation region. In this work, we will effectively eliminate the target 
fragmentation region by implementing a cut Xi;' > on the Feynman- variable representing 
the fractional longitudinal c.m.s. momentum. Target fragmentation could be accounted 
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for by transforming to the variable [29-32] 

the energies E^, Ej^ defined in the c.m.s. frame of the nucleon and the virtual photon, 
and by introducing the so-called "fracture functions" p9|. The inclusion of the latter is 
beyond the scope of this analysis and anyway not relevant numerically due to the 
cut on xp. The variable z in (^) is also better suited for dealing with corrections due to 
the finite target mass Mh- As will be demonstrated below, it is not always justified to 
neglect these. Our predictions for A production in SIDIS will therefore be made using the 
variable z. The NLO corrections to F^^^ and F^^^ in (^) can, however, be expressed 
much more conveniently in terms of zh (see Appendix C). The transformation from zh 
to z is straightforward [29-32]. 

The structure functions F^^^ and F^^^ in (^1]) are given at NLO by 

2F^'"{x,ZH,Q') = 5^e;|g(x,g^)Df(^^,g2) + ^^[g®C;®Df 

<?,<? 

+ q ® 0Df + g0 C]^ ® ] (x, ZH^ Q') | (33) 



<},<} 



+ q^C^g^n;; + g®C^g®D^ ix,ZH,Q') , (34) 



with the NLO coefficient functions Cj-^ |]20|, ^ collected in Appendix C. 



As already mentioned in the Introduction, three other possible cross sections can be 
defined when the polarization of the lepton, the initial nucleon and the hadron are taken 
into account. If both nucleon and hadron are polarized and the lepton is unpolarized, the 
expression is similar to Eqs. (pTD-(^4D above with, however, the unpolarized parton distri- 
butions and the fragmentation functions to be replaced by their polarized counterparts. 
Obviously, one also has to adapt the coefficient functions to this case: Cj^^ /S.C]-^'^^ . 
The relevant expressions can again be found in Appendix C. In the case that the lepton 
and either the nucleon or the hadron are polarized, the expression for the cross section is 
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given as in the fully inclusive case by a single structure function ' {x, zh, Q 

dAa" Aira^ 



dx dy dzfj Q"^ 
To NLO, Qi^^ can be written as 0, ^ 

2g^'"{x,ZH,Q') = $^e^|(A)g(a:,Q2)(A)Df(zH,Q 



^^"{x,Zh,Q^): 
(2-y)//^(x,^H,Q'). (35) 



+ 
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{A)q ® AC; ® (A)Df + (A)g ® A^, ® (A)Df 



+ {A)g ACi^® {A)Df]{x,ZH,Q')^ , (36) 

the position of the A and the index i = N, H depending on which particle (A^ or H) is 
polarized. Again, all NLO MS coefficient functions ACj^ are collected in Appendix C. 



Let us first turn to the entirely unpolarized case as defined in Eq. (|31|), which could 
prove invaluable for obtaining a fiavour separation of fragmentation functions not provided 
by the SIA data. Unfortunately, only three measurements of this cross section exist up to 
now m, with still rather large experimental uncertainties. It is nevertheless worth com- 
paring our predictions to the available data in order to test our proposed fragmentation 
functions in a process other than SIA. 

The original experimental results [Q are compiled in Fig. 6a, where the data are 
plotted in terms of the Feynman variable xp- As can be observed, differences between 
results at different values of the 7*p c.m.s. energy W are much larger than expected 
from the scale dependence, especially at small xp. The reason for this is simple and 
corresponds to the fact that the variable xp is not the scaling variable of this process, i.e. 
the argument of the fragmentation functions in Eq. (|3TD. As already mentioned, Eq. (|3T| ) 
should be expressed in terms of z in (|3^) pO| . At LO, z coincides with xp and also with 



the variable z^ introduced in (|3TD ||20| , |2l[| , if the mass of the A is neglected. However, it 



again turns out that, as in the case of SIA, finite-mass effects introduced by the function 



= a/1 — 4M^/ {zWy become relevant at small z for the low-VT experiments. The other 
two variables are at LO given in terms of f3 and z hj xp = (3 z and zjj = {1 + f3)/2z. 
In Fig. 6b we show the same data as in Fig. 6a, but converted to the variable z, where 
now a much better agreement between different experimental data and also with our 
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LO predictions plotted for three different typical scales can be observed. It should be 
noticed that the HI data 0] were obtained with an integrated luminosity of only 1.3 
pb~^, so a more dedicated measurement in the future will be very helpful in determining 
the unpolarized A fragmentation functions more precisely. 

The most interesting observable with respect to the determination of the polarized 
A fragmentation functions is of course the asymmetry for the production of polarized 
A's from an wnpolarized proton, defined by = g^^^ /Ff^^ with g\^^ given by 
(^) with i = H. In Fig. 7a, we show our LO and NLO predictions for HERA with 
polarized electrons and unpolarized protons using the GRV parton distributions |25| , 



integrated over the measurable range 0.1 < z < 1. The values for Q that correspond 



to each x-bin have been chosen as in [Q. Good perturbative stability of the process is 
found. As can be seen, the results obtained using the three distinct scenarios for polarized 
fragmentation functions turn out to be completely different. Since the asymmetry at small 
X is determined by the proton's sea quarks, its behaviour can be easily understood: in seen. 
1 only s quarks fragment into polarized A's, giving an asymmetry which is positive but 
about three times smaller than the one of seen. 3 where all the flavours contribute. In the 
case of seen. 2 the positive contribution from the s-quark fragmentation is cancelled by a 
negative one from u and d, resulting in an almost vanishing asymmetry. The interpretation 
is similar for the region of large x, where only the contribution involving is sizeable 
and the asymmetry asymptotically goes to J dzAD^/ J dzD^ for each scenario. 

We have included in Fig. 7a also the expected statistical errors for HERA, computed 
assuming an integrated luminosity of 500 pb^^ and a realistic value of e = 0.1 for the 



efficiency of A detection |36|. Comparing the asymmetries and the error bars in Fig. 7a 
one concludes that a measurement of at small x would allow a discrimination between 
different conceivable scenarios for polarized fragmentation functions. Fig. 7b shows our 
results vs. z for fixed x = 5.6 ■ 10^^. Again, very different asymmetries are found for 
the three scenarios. In this plot we also include the expectation from scenario 1 with a 
maximal gluon polarization (see Sec. 3) which, as expected, predicts a larger asymmetry 
at small z (comparable to the one of scenario 3) mainly due to the contribution of the 
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radiated sea. 



In Fig. 8 we show the same observable for the case of HERMES, where the values 
were chosen as for the inclusive DIS measurements by HERMES |37|. This fixed target 



experiment analyses a different kinematical region of both larger x and z i^z > 0.3) and 
hence could provide complementary information. As can be seen by comparing Figs. 7b 
and 8b, the asymmetry for scenario 3 shows a similar behaviour for both experiments, 
which is expected as all the fragmentation functions are equal and the parton distributions 
cancel in the ratio (at LO). However, the asymmetries predicted by scenarios 1 and 2 
change quite a bit when going from HERA collider to fixed target energies due to the 
fact that for the values of x probed at HERMES the contributions from the valence 
distributions dominate. Thus more weight is given to the fragmentation functions, 
and the contributions involving are suppressed, in contrast to the situation for the 
small x-region to be explored by HERA. Again we show in Fig. 8b also the results for 
scenario 1 with a maximally saturated gluon fragmentation function at the input scale 
which leads to results hardly different from the standard one since the sea contribution is 
negligible at large z. 

Finally, the particular case of both target and hadron being polarized was originally 
proposed as a very good way to obtain the As distribution [R3] . The underlying assump- 



tion here was that only the fragmentation function AZ)^ is sizeable (as realized, e.g. in 
our scenario 1), and that therefore the only contribution to the polarized cross section has 
to be proportional to AsAD^. In order to analyze the sensitivity of the corresponding 
asymmetry to As, we compute it using the two different GRSV sets of polarized parton 



densities of the proton |26|, which mainly differ in the strange distribution: the so-called 
"standard" set assumes an unbroken SU{3)f symmetric sea, whereas in the "valence" 
scenario the sea is maximally broken and the resulting strange quark density is quite 
small. 

The results for HERA are shown in Fig. 9. Unfortunately - and not unexpectedly 
- it turns out that the differences in the asymmetry resulting from our different models 
for polarized A fragmentation are far larger than the ones due to employing different 
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polarized proton strange densities. In addition, a distinction between different As would 
remain elusive even if the spin-dependent A fragmentation functions were known to good 
accuracy, as can be seen from the error bars in Fig. 9 which were obtained using the same 
parameters as before. 

5 Summary and Conclusions 

We have performed a detailed QCD analysis of the production of A baryons in e'^e~ 
annihilation and semi-inclusive deep-inelastic scattering. 

Working within the framework of the radiative parton model, our starting point has 
been a fit to unpolarized data for A production taken in e^e~ annihilation, yielding 
a set of realistic unpolarized fragmentation functions for the A. We have then made 
simple assumptions for the relation between the spin-dependent and the unpolarized A 
fragmentation functions at the input scale for the Q^-evolution. Taking into account the 
sparse LEP data on the polarization of A's produced on the Z-resonance, we were able to 
set up three distinct "toy scenarios" for the spin-dcpcndcnt A fragmentation functions, to 
be used for predictions for future experiments. We emphasize that our proposed sets can 
by no means cover all the allowed possibilities for the polarized fragmentation functions, 
the main reason being that the LEP data are only sensitive to the valence part of the 
polarized fragmentation functions. Thus, there are still big uncertainties related to the 
"unfavoured" quark and gluon fragmentation functions, making further measurements in 
other processes indispensable. 

Under these premises, we have studied A production in semi-inclusive dccp-inelastic 
scattering. Existing data for the production of unpolarized A's are well described by our 
fragmentation functions determined from the e"'"e~ annihilation data. Turning to spin 
transfer asymmetries sensitive to the longitudinal polarization of the produced A's, we 

— * — * 

have considered both ep — > AX and ep — > AX scattering. It turns out that in the first case 
SIDIS measurements at HERA (with spin- rotators in front of the HI and ZEUS detectors) 
and at HERMES should be particularly well suited to yield further information on the 
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ADj : differences between the asymmetries obtained when using different sets of are 
usually larger than the expected statistical errors. In contrast to this, having a polarized 
proton target (or beam) does not appear beneficial as far as A production is concerned. 

A Fortran package containing our unpolarized and polarized LO and NLO A frag- 
mentation functions can be obtained by electronic mail from Daniel. Deflorian@cern.ch, 
Marco.Stratmann@durham.ac.uk, or Werner. Vogelsang@cern.ch upon request. 
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Appendix A: Unpolarized SIA 



The NLO (MS) coefficients C]^ in (0) and (|) are given by 0, |21[: 

'ln(l-2)\ 3 1 



+2 



{l + z' 



1 - z 



2(1-^)+ 



Cl{z) 
Ciiz) 
Ciiz) 



2Cp 
4Cf 



_^ln. + ^(l-.)+Q.^J).(l-.) 
l±il^ln(.^(l-.))-2l^ 



(A.l) 



(A.2) 
(A.3) 
(A.4) 



with Cf = 4/3. Note that in the expressions for C^g we have taken the factorization 
scales for the final-state mass singularities to be equal to the hard scale Q of the process, 
as we did in all our numerical applications. The "+" -prescription is defined as usual by 



\zf{z){g{z))^^ f dz[f{z)-f{l)]g{z) . 

^0 



(A.5) 



The electroweak charges in 



are given by 



el = el- 2e,XimVeV, + X2(Q')(1 + + 



(A.6) 
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where 



X2[S) 



s{s-Ml) 



16 sin^ GvF cos2 {s - M^f + r|M| 
1 



(A.7) 



256 sin^ Qw cos^ {s - M^f + r|M| ' 
Here Cq is the fractional electromagnetic quark charge, and Mz and T z are the mass and 
the decay width of the Z boson, respectively. The other electroweak couplings are given 
in terms of the Weinberg angle Qy/ by 



Vu 



-1 + 4sin^ Qw 
g 

+1 - - sin^ Qw 

-1 + - sm (dw 

o 



(A. 



Appendix B: Polarized SI A 



The NLO MS coefficients ^0]^^ in (0)-(^ read: 



ACliz) 



<? 

2Cf 



{2-z) \n{z^{l-z)) -4 + 3, 



AC'(z) = Cliz) 



where the effective charges on the Z-resonance are given by 



9a 



2x2(M|)K(l + K') 



(B.l) 
(B.2) 
(B.3) 
(B.4) 

(B.5) 
(B.6) 



where Au = —Ad = 1 and Ve, Vq, X2{s) have already been defined in ( |A.7| ), ( [A.8| ). The 
structure functions and in Eqs. (0), ( pQ]) are purely non-singlet and therefore do 
not receive a gluonic correction. 



One should note that the NLO quark corrections for the unpolarized case, see Eqs. 
( |A.1| ) and ( [A.3[ ), and the ones for the spin-dependent parity violating structure functions 
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gf^ and in ( |B.3| ), ( [B.4| ) are identical, which resuhs from identical tensorial structures 
at the parton level. The expressions for /\Cq{z) and /\C]^{z) in the MS scheme were 
already derived in |l^]. The difference /S.C]^{z) — C^{z) in ( p.l|) is independent of 
the regularization prescription chosen and coincides with the one found in pSf by using 
off-shell gluons to regularize the coUinear singularities. 



Appendix C: Unpolarized and Polarized SIDIS Coef- 
ficient Functions 



Here we list all unpolarized and polarized NLO (MS) coefficients (A)C.'.'.' for SIDIS as 
introduced in Section 4. To keep the expressions as short as possible it is convenient to 
define the following abbreviations 



^i(0 = (i+e 



P,M) = + {-^ - if 
2^ [H^-iY 



if 



AP,,{o = e - {I - if 



2-e , 



V 1-e 



^2(0 



Ine . 



(C.l) 



Note that in what follows we always suppress the argument {x, z) of the coefficient func- 
tions. M and Mp denote the factorization scales for initial and final state mass singular- 
ities, respectively. Note that for all our numerical calculations we have chosen as usual 
M = Mp = Q- AH results presented here are given in the MS scheme, and in case of 
the spin-dependent coefficients AC*g the additional finite subtractions that are required 



when using the 75 prescription of have been performed along the lines discussed in 



Coefficients for eA^ ^ e'HX: 



^1 

^qq 



Cf - 8(5(1 + 

5{l - x) \p,,{z) In ^ + L,{z) + L^{z) + {I - z) 



+ 
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^gq 



C. 



C. 



6{l-z) 
1 



+ 



1 + z 



1 + x 



{l-x)+{l-z)+ (1-^)+ 

1 + x 



+ 2{l + xz) 
+ 



zS{l -x) + 2{l + x- xz) - 



^qg 



8{\ - z) 



Pgg{x) In 



1-x 

M2 



X 



qg 



Pqgi^j ' (1 - ■ z 
ACpxz 

4:Cfx{1 - z) 

4x(l - x) 



+ --2 



+ 2x(l -x) 



+ 



Coefficients for eTV ^ e'HX: 





= Ci,-2CH1 






= C^i^-2C^z(l-a;) 






^) 


AP,,(x)ln(^ 








1 1 




[(1 


^ + -- 



M2 X 



+ 2(l-x) 



+ 



Coefficients for eTV — > e'HX: 



AC, 



99 



Aa^ = 



99 



AC^^^ + 2Cf (1-^)5(1-3;) 
Ci.|AP,,(^) 5(1 -x) In 



01 



^(1-^) + 



2(1 - ^)5(1 - x) - 2(1 + X - ^) + 



1 + a; 



^^qg ~ ^qg ~ ^qgi^) 



1-z 



Coefficients for eN — > e'HX: 



AC. 



1,NH 

qq 



Clg + 2CF{l-z)d{l-x) 
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AClf"" = AC^^-2CMl'^) (C.15) 

AClf^ = AC^^ - AP,,{x)^—^ (C.16) 

AC^f" = (C.17) 

AC^gf^ = -4Cfx{1-z) (C.18) 

ACi;-^^ = . (C.19) 



We note that all our results for the spin-dependent coefficients in (C.8)-(C.19) co- 



incide with the ones presented in and also fully agree with the results of |3T|, ^ 
if one carefully disentangles in these papers the contributions from the current and the 
target fragmentation regions (in the same way our unpolarized results in (C.2)-(C.7) are 



in agreement with ref. pO[). In addition, one has to account for the slightly different 



factorization scheme used in ref. 31 



Finally let us show how to deal with the - distributions appearing in the expres- 
sions above. The - distribution was already defined in Eq. ( |A.5D in Appendix A. In 
practice, however, the lower limit of the integration in ( |A.5| ) is different from zero, hence 



the distributions have to be modified according to 

1 1 



(1-0 
Inl-e 



Inl-e 



+ ln(l-A)5(l-0, 



(C.20) 



where ( )yi is defined as in ( [A.5| ) but with the lower integration limit replaced by A. In 
addition, in the coefficients (A)C*^ listed above also double - distributions appear, 
which can be defined in analogy with Eq. ( |A.5|) by 



dx / dz- 



X] 



[l-z) 



' f\,^ f{x,z)-f{l,z)-f{xA) + f{l,l) 
Jo 



(C.21) 



Again, in practice the lower integration limits are both different from zero, say, A for 
the X integration and B for the z integration, and the distribution defined above can be 
rewritten as 

^ ^ ^ \n{l - B)6{1 - z) + 



(l-x)+(l-^)+ {1~x)a{1-z)b (1-x), 
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— ln(l - A)d{l -x)+ ln(l - A) ln(l - B)6{1 - x)6{l - z) . (C.22) 

(1 - z)b 
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Figure Captions 

Fig. 1 Comparison of our LO and NLO results for (1/ atot)d(y / dxE according to Eq. (^) 
with all available data on unpolarized A production in e^e~ annihilation [^ ^. Note 
that only data points with xe > 0.1 have been included in our fit (see text). 

Fig. 2 z-dependence of our LO and NLO fragmentation functions as specified in Eq. (|23|) 
and Table 1 at = iqo and Q"^ = 10^ GeV^. 

Fig. 3 Q^-dependence of our LO and NLO q and g fragmentation functions at fixed 
z = 0.05, 0.1, 0.5. 

Fig. 4 Comparison of LEP data [10-12] and our LO and NLO results for the asymmetry 
in (|26|), using the three different scenarios as described in the text. 

Fig. 5 LO and NLO partonic fragmentation asymmetries Aj = AD^/D^ for f = s, 
u = d, and g at = 10 GeV^. In the LO plot we also show (dot-dashed lines) the 
effect of assuming a maximally polarized gluon distribution at the initial scale for 
scenario 1. The NLO results in this case are very similar to the LO ones and are 
therefore not shown. 

Fig. 6 a) Compilation of the original SIDIS data [^] in terms of xp- b) Comparison of 
our LO predictions with the data converted to the variable z (see text). 

Fig. 7 LO and NLO predictions for the SIDIS asymmetry for unpolarized protons and 
polarized A's and leptons (see text) for our three distinct scenarios of polarized 
fragmentation functions. In a) we also show the expected statistical errors for such 



28 



a measurement at HERA, assuming a luminosity of 500 pb~^, a beam polarization 
of 70%, and a A detection efficiency of 0.1. In b) we include the expectation for 
scenario 1 with a maximal gluon polarization at the initial scale. 

Fig. 8 The same as in Fig. 7, but now for HERMES kinematics. 

Fig. 9 LO and NLO predictions for the SIDIS asymmetry for polarized protons but 
unpolarized leptons for two different sets of polarized parton distributions taken 
from pGf . Also shown are the expected statistical errors for such a measurement at 



HERA, calculated for the parameters already used for Fig. 7. 
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